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Abstract. We report precise Doppler measurements of HD 142022 obtained during the past six years with the CORALIE 
echelle spectrograph at La Silla Observatory together with a few additional observations made recently with the HARPS echelle 
spectrograph. Our radial velocities reveal evidence of a planetary companion with an orbital period P = 1928^3, days, an 
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and a velocity semiamplitude K = 92+29 m s . The inferred companion minimum mass is M2 sin / : 



Mjup and the semimajor axis a = 3.03 ± 0.05 AU. Only one full orbital revolution has been monitored yet, and the 
periastron passage could not be observed since the star was too low on the horizon. The eccentricity and velocity semiamplitude 
remain therefore quite uncertain and the orbital solution is preliminary. HD 142022 is a chromospherically inactive KO dwarf, 
metal rich relative to the Sun, and is the primary component of a wide binary. HD 142022 b is thus a new "planet in binary" 
candidate, and its high eccentricity might be due to secular interactions with the distant stellar companion. 
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The CORALIE survey for southern extrasolar planets has been 
going on since June 1998. This high-precision radial-velocity 
programme makes use of the CORALIE fiber-fed echelle spec- 
trograph mounted on the L2-m Euler Swiss telescope at La 
Silla Observatory (ESO, Chile). The sample of stars monitored 
for the extrasolar planet search is made of 1650 nearby G and 
K dwarfs selected according to distance in order to have a 
well-defined volume-limited set of stars cUdry et al. 2000). The 
CORALIE sample can thus be used to address various aspects 
of the statistics of extrasolar planets. 

Stellar spectra taken with CORALIE are reduced online. 
Radial velocities are computed by cross-correlating the mea- 
sured stellar spectra with a numerical mask, whose nonzero 
zones correspond to the theoretical positions and widths of 
stellar absorption lines at zero velocity. The resulting cross- 
correlation function (CCF) therefore represents a flux weighted 
"mean" profile of the stellar absorption lines transmitted by the 
mask. The radial velocity of the star corresponds to the mini- 
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mum of the CCF, which is determined by fitting the CCF with 
a Gaussian function. The initial long-term velocity p recision of 
CORALIE was about 7 ms ' dOueloz et alJl200ol) . but since 
2001 the instrumental accuracy combined with the simultane- 
ous Th Ar-reference technique is better than 3ms' (Pepe et al.j 
l2002al) . This implies that for many targets the precision is now 
limited by photon noise or stellar jitter. 

After seven years of activity, the CORALIE survey has 
proven to be very successful with the detection of a signifi- 
cant fraction of the known extrasolar planetary candidates (see 
the previous papers in this series). As the survey duration in- 
creases, new planetary candidates with orbital periods of sev- 
eral years can be unveiled and their orbit characterized. This is 
the case of the companion orbiting HD 142022: with a period 
longer than 5 years it has just completed one orbital revolu- 
tion since our first observation in 1999. This growing period- 
interval coverage is very important with regard to formation 
and migration models since observational constraints are still 
very weak for periods of a few years. 

HARPS is ESO's High-Accuracy Radial-Velocity Planet 
Searcher ( Pepe et al...2002b. ,2004; Mayor et al...2003.) . a fiber- 
fed high-resolution echelle spectrograph mounted on the 3.6- 
m telescope at La Silla Observatory (Chile). The efficiency 
and extraordinary instrumental stability of HARPS combined 
with a powerful data reduction pipeline provides us with very 
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Table 1. Observed and inferred stellar parameters for 
HD 142022 (see text for references). 



Parameter 


Unit 


Value 


Spectral Type 




KO 


V 


(mag) 


7.70 


B-V 


(mag) 


0.790 


n 


(mas) 


27.88±0.68 


My 


(mag) 


4.93 




(K) 


5499±27 


log 5 


(cgs) 


4.36±0.04 


[Fe/H] 


(dex) 


0.19±0.04 


L 


(Lo) 


1.01 




(Mo) 


0.99 


V sin ; 


(kms-') 


1.20 


log(^HK) 




-4.97 



high precision radial-velocity measurements, allowing the de- 
tectio n of companions of a few Earth masses around solar-type 
stars (^S antos et al.ll2004l) . Benefiting from this unprecedented 
precision, a part of the HARPS Consortium Guaranteed-Time- 
Observations programme is devoted to the study of extraso- 
lar planets in a continuation of the CORALIE survey, to al- 
low a better characterization of long-period planets and multi- 
ple planetary systems. HD 142022 is part of this programme, 
and the few HARPS measurements obtained so far already con- 
tribute to improve the orbital solution based on CORALIE data. 

The stellar properties of HD 142022 are summarized 
in Sect. |5] Section |3 presents our radial-velocity data for 
HD 142022 and the inferred orbital solution of its newly de- 
tected companion. These results are discussed in Sect.|3 show- 
ing that the planetary interpretation is the best explanation for 
the observed velocity variation. Our conclusions are drawn in 
Sect.|5] 

2. Stellar characteristics of HD 142022 

HD 142022 (HIP 79242, Gl 606.1 A) is a bright KO dwarf in 
the Octans constellation. The astrometric parallax from the 
Hipparcos catalogue, n = 27.88 ± 0.68 mas (ESA 1997), sets 
the star at a distance of 36 pc from the Sun. With an appar- 
ent magnitude V - 7.70 (ESA 1997) this implies an absolute 
magnitude of My - 4.93. According to the Hipparcos cata- 
logue the color index for HD 142022 is fl- V = 0.7 90. Using a 

' ■ 1 

bolometric correction BC - -0.192 (Flower 1996) and the so- 
lar absolute magnitude Mq°' - 4.746 (Leieune et al...l998) we 
thus obtain a luminosity L - 1.01 Lq. The stellar parameters 
for HD 142022 are summarized in Tabled 

A detailed spectroscopic analysis of HD 142022 was per- 
formed using our HARPS s pectra in order to o btain accurate 
atmospheric parameters (see lSantos et al ] 42005) for further de- 
tails). This gave the following values: an effective temperature 
Teff = 5499 + 27 K, a surface gravity log g = 4.36 + 0.04, and a 
metallicity [Fe/H] - 0.19 ± 0.04. Using these param eters and 
the Geneva stellar evolution code jMevnet & Maederi.200 0) we 
deduce a mass M* - 0.99 Mq. According to evolutionary mod- 




I ' — ' ' ' ' ' ' — ' — ' ■ ' — ' ' ' ' — ' ' ' ' — I 

3964 3966 3968 3970 3972 3974 

Wavelength [A] 



Fig.l. Can H (/I = 3968.47 A) absorption Une region of the 
summed HARPS spectra for HD 142022. 



els, HD 142022 is an old main-sequence star, in agreement with 
the KOV spectral type quoted in the Hipparcos catalogue. 

The cross-correlation function can be used to derive stel- 
lar quantities affecting line profiles such as the projected rota- 
tional velocity. From the CORALIE spectra we derive v sin / - 
1.20 kms"' dSantos et al...2002) . Combining this result to the 
stellar radius given by the best evolutionary model {R - 
1.15 Rq) we obtain an upper limit of 48 days for the rotational 
period. 

From the HARPS spectra we can compute the log(^^j^) 
activity index measuring the chromospheric emission flux in 
the Can H and K lines. This index is a useful estimator of the 
radial-velocity jitter that can be expected from intrinsic stellar 
variability. Figure ^ shows the Ca ii H absorption line region 
for HD 142022. No emission peak is visible at the center of the 
absorption line, indicating a rather low chromospheric activity. 
This is corroborated by the log(/?^j,) value of -4.97, typical of 
inactive stars. 

HD 142022 is part of a wide binary. Its companion, 
LTT 6384 (Gl 606. IB), is a late-K star about 22 arcseconds 
away and with an apparent magnitude V - 11.2. The two 
stars are listed in the NLTT and LDS catalogs (lT,iivtenlll94nl 
[1979) indicating very similar proper motions. They were also 
observed by the Hipparcos satellite, but the proper motion 
of LTT 6384 could not be determined. The apparent sep- 
aration of the pair, nevertheless, remained close to 22 arc- 
seconds from 1920 to 2000 which, given the proper mo- 
tion of HD 142022 (fi,, cos 6 = -337.59 + 0.60 masyr"', 
fxg - -31.15 + 0.72 masyr"'), is an indication that the pair 
is indeed a bound system. This conclusion is strengthened 
by the fact that the CORAVEL radial velocitie s of the tw o 
stars are identical within uncertainties ( Nordstro m et al.l2004l) . 
Using the positions given by the Tycho-2 catalogue and its 
supplement-2 (ESA 1997), we obtain a projected binary sep- 
aration of 820 AU. This translates into an estimated binary 
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Fig. 2. Lomb-Scargle periodogram of the CORALIE radial ve- 
locities for HD 142022 (top) and the corresponding window 
function (bottom). The power is a measure of the statistical sig- 
nificance of the signal, not of its true amplitude. 

semimaior axis of 1033 AU, using the relation ajr - 1.26' 
jFischer & Marc\ll992l) . 

In the Hipparcos catalogue, HD 142022 is classified as an 
unsolved variable and is suspected of being a non-single star. 
Indeed, a Lomb-Scargle periodogram of the Hipparcos pho- 
tometry shows no clear signal standing out, but some extra 
power spread over many frequencies, especially at short pe- 
riod (few days). Performing a detailed study of the Hipparcos 
photometry for HD 142022 is beyond the scope of this paper, 
and is not fundamental since the periods involved are much 
shorter than that of the signal detected in radial velocity. We 
will nonetheless briefly come back to this issue in Sect. |3 be- 
cause the potential non-single status of the star may be a con- 
cern. 



3. Radial velocities and orbital solution 

HD 142022 has been observed with CORALIE at La Silla 
Obervatory since July 1999. Altogether, 70 radial-velocity 
measurements with a typical signal-to-noise ratio of 25 (per 
pixel at 550 nm) and a mean measurement uncertainty (in- 
cluding photon noise and calibration errors) of 9.7 ms ' were 
gathered. HD 142022 is als o part of the H ARPS high-precision 
radial-velocity programme (iLovis et a l. 2005) and, as such, was 
observed 6 times between November 2004 and May 2005. 
These observations have a typical signal-to-noise ratio of 90 
and a mean measurement uncertainty of 1 m s ' . 



' Strictly speaking, the relation used here to translate the projected 
separation of a wide binary into a semimajor axis is valid only statis- 
tically. It can thus be highly inaccurate for an individual system. 



Fig. 3. CORALIE and HARPS radial velocities for HD 142022. 
CORALIE velocities are shown as dots, HARPS data are plot- 
ted as circles. As each instrument has its own zero point, a 
velocity off'set of -8.8 ms ' has been added to the HARPS 
data (see text for further details). Error bars for HARPS data 
(1 m s"') are the same size as the symbols. 

The root-mean-square (rms) of the CORALIE radial ve- 
locities is 26.3 ms"', indicating some type of variability. The 
Lomb-Scargle periodogram of these velocities is shown in 
Figure 121 The highest peak coiTesponds to a period of 1926 
days, which is clearly visible in the plot of our radial veloci- 
ti es as a function of time (Fig.|5}- Using the expressions given 
in lScargl3 (Il982 !). the false alarm probability for this signal is 
close to 10"^. This low value was confirmed using Monte Carlo 
simulations, in which data sets of noise only were generated 
with velocities drawn at random from the residuals around the 
mean. None of the 10^ simulated data set exhibited a maximum 
periodogram power exceeding the observed value, yielding a 
false alarm probability <10"^. Figure |4] shows the CORALIE 
and HARPS radial velocities phased with a similar period and 
the coiTesponding best-fit Keplerian model. The resulting or- 
bital parameters are P = 1928 days, e - 0.53, K - 92 ms"', 
implying a minimum mass M2sin/ = 5.1 Mjup orbiting with 
a semimajor axis a - 3.03 AU. The orbital elements for 
HD 142022 are listed in Table|5] Since each data set has its own 
velocity zero point, the velocity off'set between the two instru- 
ments (HARPS and CORALIE) is an additional free parameter 
in the fitting process. Note that the first two HARPS measure- 
ments are contemporary with the last CORALIE observations 
(Fig-E), and the curvature seen in the HARPS data fits perfectly 
well the CORALIE velocities taken during the previous revo- 
lution (Fig.0}. The velocity offset is therefore well constrained 
though the number of HARPS measurements is small. 

As can be seen in Fig.|3 our data span only one orbital pe- 
riod and the phase coverage is very poor near periastron since 
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the star was too low on the horizon to be observed at that time. 
This is why the orbital eccentricity and the velocity semiampli- 
tude are poorly constrained. The orbital solution is thus pre- 
liminary and additional measurements taken during the next 
revolutions will be needed to obtain more accurate orbital pa- 
rameters. It should be noted that for such an eccentric orbit the 
semiamplitude K is strongly correlated with the eccentricity. 
As our present analysis is more likely to have overestimated the 
eccentricity, the semiamplitude K might in fact be smaller, im- 
plying a smaller minimum mass for the companion. The com- 
panion minimum mass is thus very likely to be in the planetary 
regime, whatever the exact eccentricity. 

The uncertainties in the orbital parameters were determined 
by applying the fitting technique repeatedly to many sets of 
simulated data. 1000 simulated data sets were thus constructed 
by adding a residual value (drawn at random from the resid- 
uals) to the best-fit velocity corresponding to each observing 
time. For each realization, the best-fit Keplerian orbit was de- 
termined. This has been done by using the period obtained 
from the Lomb-Scargle periodogram as an initial guess for the 
Keplerian fit. We then used a local minimization algorithm to 
find the best fit, trying several initial starting values for T, u 
and K. The quoted uncertainties correspond to the Icr confi- 
dence interval of the resulting set of values for each orbital pa- 
rameter. 

The rms to the Keplerian fit is 10.8 m s ' for the CORALIE 
data and 1 .4 m s ' for the HARPS measurements, yielding a re- 
duced of 1.5. The Keplerian model thus adequately explains 
the radial-velocity variation, though for both instruments the 
rms is slightly larger than the mean measurement uncertainties. 
A periodogram of the velocity residuals after subtracting off 
the best-fit Keplerian shows no signal with significant power 
anymore. Furthermore, no changes in stellar line profiles (as 
quantified by the CCF bisector span, see Oueloz et al. (2001) ) 
are seen in our data. 

4. Discussion 

So far, it has been implicitly assumed that the radial-velocity 
variation observed for HD 142022 stems from the presence of 
a planetary companion. But it is well known that other phe- 
nomena can induce a periodic variation of the radial veloc- 
ity of a star. Activity related phenomena such as spots, plages 
or inhomogeneous convection are such candidates. They can 
however not be put forward to explain the signal observed for 
HD 142022 since the period is 1928 days, much too long to be 
related in any way to the rotational period of the star Stellar 
pulsations may also cause radial-velocity variations, but no 
known mechanism could be invoked to sustain large-amplitude 
long-period oscillations in a KO main-sequence star. Stellar 
pulsations can therefore be ruled out as well. 

The presence of a second and varying faint stellar spec- 
trum superimposed on the target spectrum can also induce spu- 
rious radial- velocity variations mimicking a pla netary-ty pe sig- 
nature (see the semi nal case of HD 41004, Sant os et al.l 112002): 

^^^^^^^^^^^^^^^^^^^^ ^^^^^^^^1 ^j^^H^H^HHH^ 

IZucker et alJ (120031) '). In such a case, the target is the unrecog- 
nized component of a multiple stellar system. Given the fact 
that HD 142022 was suspected of being non-single on the ba- 
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Fig. 4. Phased CORALIE and HARPS radial velocities for 
HD 142022 (top). CORALIE observations are shown as dots, 
HARPS data are plotted as circles. The solid line is the best- 
fit orbital solution. Residuals to the Keplerian fit are also dis- 
played (bottom). Error bars for HARPS data are the same size 
as the symbols. 

sis of Hipparcos data, this is a possibility to take seriously into 
consideration. The evolution of the cross-correlation bisector 
span as a function of radial velocity can be a useful tool in dis- 
entangling planetary signatures fr om spu rious radial-velocity 
signals (Oueloz etal. 2001; Santo s et'aD Eo02). For a signal 
of planetary origin the bisector span is constant whatever the 
value of the radial velocity, whereas for blended systems it is 
often correlated with the mea s ured r adial velocity (see for ex- 
ample Fig. 3 of ISantos et al.l ( 120021) '). No such correlation is 
visible for HD 142022. We have also searched for the pres- 
ence of a second spectrum in our CORALIE spectra using 
multi-orde r TODCOR, a two-dime nsional cross-corre lation al- 
gorithm ( Mazeh & Zucker..l994.: iZucker et alJl2003l) . but did 
not find anything convincing. These negative results do not al- 
low us to formally discard the blend scenario, but they render 
it an unlikely possibility. To sum up, the 1928-day signal ob- 
served in our radial velocities for HD 142022 is most likely to 
be caused by the gravitational perturbation of a planet orbiting 
the star. 

The Keplerian fit to the radial velocities of HD 142022 
implies that the stellar orbit has a semimajor axis cii sin; - 
0.015 AU. Given the stellar parallax, this translates into an 
angular semimajor axis a\ sin; — 0.41 mas. If the true mass 
of HD 142022 b were much larger than its minimum mass of 
5.1 Mjup, the stellar wobble might be detected in the Hipparcos 
astrometry. This potential wobble may, however, be partially 
absorbed into the solution for proper motion and parallax since 
the orbital period is longer than the 2.7-year duration of the 
Hipparcos measurements. We searched for an astrometric wob- 
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Table 2. Orbital parameters for HD 142022 b. 



Parameter 


Unit 


Value 


P 


(days) 


1928!^^ 


T 
e 
7 

(0 


(JD-2400000) 

(m s"') 
(deg) 


50941+^0 

"•-'-'-0.18 

q 7qo+0.007 
170+* 


K 


(m s-') 


Q9+102 


M2 sin i 
a 


(Mj„p) 
(AU) 


5 r- * 


Velocity offset (HARPS) 


(m s-') 


-8.8!?i 


A'nieas (CORALIE+HARPS) 




70+6 


rms (CORALIE) 


(m s-') 


10.8 


rms (HARPS) 


(m s-') 


1.4 



ble in the Hipparcos data for HD 142022, but did not find any- 
thing significant. 

HD 142022 b is one of the planets with the longest orbital 
period found in a wide binary so far, but its separation of 3 AU 
is still very small compared to the estimated binary semimajor 
axis of 1033 AU. HD 142022 b thus orbits well inside the sta- 
bility zone, whatever the e xact orbital parameters of the binary 
llHoknan & Wiegertlll999l) . The presence of a distant stellar 
companion may nonetheless cause significant secular pertur- 
bations to the planetary orbit. In particular, if the planetary or- 
bital plane is inclined relative to the binary plane, the planet can 
undergo large-amplitude eccentricity oscillations due to t he so- 
called Kozai mechanism (Kozai (1962); see also Holman et af] 
lll997l) : llnnanen et alJ ( ll997l) : lMazeh et aP (Il997h '). The Kozai 
mechanism is effective at very long range, but its oscillations 
may be suppressed by other competing sources of orbital per- 
turbations, such as general relativity effects or perturbations re- 
sulting from the presence of an additional companion in the 
system. Regarding HD 142022, we ha ve estimated the ratio 
^'Kozai/^'cR using equations 3 and 4 of iHolman et alJ lll997h 
with the values - 1 / V2 and Ms - 0.6 M© for the binary 
eccentricity and secondary component mass, respectively. This 
yields f Kozai = 1-25 10^ years and PKozai/^'cR = 0.35, indicat- 
ing that Kozai oscillations could take place is this system, since 
their period is shorter than the apsidal period due to relativis- 
tic effects. Although not well constrained, the eccentricity of 
HD 142022 is clearly quite high, and such a high eccentricity 
is not surprising if the system undergoes Kozai oscillations. 

5. Conclusion 

We report a 1928-day radial-velocity variation of the KO dwarf 
HD 142022 with a velocity semiamplitude of 92 ms"'. From 
the absence of correlation between stellar activity indicators 
and radial velocities, and from the lack of significant spectral 
line asymmetry variations, the presence of a planetary compan- 
ion on a Keplerian orbit best explains our data. The Keplerian 
solution results in a M2sini - 5.1 Mjup companion orbiting 
HD 142022 with a semimajor axis a - 3.03 AU and an ec- 
centricity e - 0.53. Although HD 142022 b orbits the pri- 
mary component of a wide binary, its characteristics, includ- 



ing minimum mass and orbital eccentricity, are typical of the 
long-period planets found so far around G and K dwarfs. 

One of the most surprising properties of extrasolar plan- 
ets revealed by ongoing radial-velocity surveys is their high 
orbital eccentricities, which challenge our current theoretical 
paradigm for planet formation. Several mechanisms have thus 
been proposed to account for eccentric planetary orbits. One of 
them is the Kozai mechanism, a secular interaction between a 
planet and a wide binary companion in a hierarchical triple sys- 
tem with high relative inclination. Although the Kozai mech- 
anism can be put forward to explain the high eccentricity of 
a few planetary companions: 16 Cyg Bb ("Holman etalJ ll997l: 
^azeh et al. 1997), HD 80606 b ( Wu & Murray 2003) and pos- 
sibly HD 142022 b, it seems impossible to explain the observed 
eccentricity distribution of extrasolar planets solely by invok- 
ing the presence of binary companion s Ijakeda & Rasio 200^. 
According to lTakeda & Rasiol ( l2005l) . Kozai-type perturbations 
could nonetheless play an important role in shaping the eccen- 
tricity distribution of extrasolar planets, especially at the high 
end. In this regard, ongoing programmes aiming at searching 
for new (faint) companions to stars with known planetary sys- 
tems, or aiming at estimating the frequency of planets in binary 
systems should soon bring new observational material, and en- 
able us to refine our present knowledge. 
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